Pulsed light beams in vacuum with superluminal and negative group velocities 
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Gouy's phase of transversally limited pulses can create a strong anomalous dispersion in vacuum 
leading to highly superluminal and negative group velocities. As a consequence, a focusing pulse 
can diverge beyond the focus before converging into it. A simple experiment is proposed. 



Propagation of light pulses at superluminal velocities 
has received a good deal of attention in recent years 
|fl, 0, [|. Superluminal group velocities in light pulses 
have been predicted and experimentally demonstrated in 
evanescent modes of undersized waveguides [[l], in X- 
Besssel waves ||, and in pulses travelling in transpar- 
ent materials with pronounced enough anomalous disper- 
sion H , one of the most striking results being a negative 
superluminal group velocity of an undcformable pulse, 
which in practice means that the pulse exits the disper- 
sive material before entering ||. 

In this letter we show that arbitrarily high and negative 
group velocities in temporally undeformable pulses are 
also possible in vacuum as an effect of the phase anomaly 
j|, or Gouy's phase associated to diffraction. 

In fact, slight superluminal behavior in collimated |j| 
and focused [[| pulsed Gaussian beam has been recently 
described. Though these results were obtained under the 
approximate paraxial theory of light beam propagation 
, they have been shown to be correct and accurate from 
the nonparaxial vectorial Kirchhoff-Sommerfeld diffrac- 
tion formula || ||. 

Discrepancy of the group velocity from c in arbitrary 
pulsed light beams is shown here to be a consequence of 
the frequency dependence of Gouy's phase, a dispersion 
which can be controlled by properly selecting the pulse 
transversal profile. We then find it rather simple, by us- 
ing standard optics such as a lens and a radially graded 
mirror ]9|, [To), to generate a strong enough anomalous 
dispersion so as to produce arbitrarily high superluminal 
and negative group velocities in a focusing pulse, whose 
temporal form is moreover found to remain nearly un- 
changed. 

Let us start by considering the three-dimensional wave 
packet, or pulsed light beam 

1 f°° 

E(r,z,t) = — / dui Eu(r, z) exp(— iut), (1) 
it Jo 

r = (x 2 + y 2 ) 1 / 2 , that results from the superposition of 
monochromatic light beams E u {r, z) of different frequen- 
cies ui. 



All components E u (r,z) are assumed to be parax- 
ial about the positive z-axis, and further to have, for 
simplicity, revolution symmetry around this axis. The 
slowly varying complex amplitude V>a>(Vj z) of the parax- 
ial light beam E u {r,z) [i.e., E u — ifj u exp(iwz/c), with 
c the vacuum light speed of plane waves] obeys the 
paraxial wave equation Aj_-0 w + 2i(u> / c)d z ij) u = 0, where 
Aj_ = d xx + d yy = (l/r)d r (rd r ) is the transversal Lapla- 
cian. Writing ip u = exp (i^ a ,), i.e., in terms of its 
real amplitude a u {r, z) = y/ I u (r, z) > [squared root 
of intensity I u {r, z)\ and phase ^(r, z), it is a straigh- 
forward calculation to obtain the propagation equa- 
tion for the on-axis (r = 0) phase 4> w {z) = 4> u (0,z), 
or Gouy's phase, as d z <j> u {z) — (c/2oj)C Ld (z), where 
Cuj(z) = Aj_a w (r, z)| r =o/ a w(0, z), and we have assumed 
that the light beam is smooth enough at on-axis points 
[d r (j)(r, z)\ r= o = 0]. The axial variation of the Gouy's 
phase d z (f>(z) of a monochromatic light beam can thus 
be inferred from the property of the transversal inten- 
sity profile C LJ (z), an easily measurable quantity having 
an intuitive meaning: If the normalized intensity pro- 
file is expanded in power series about r — 0, we obtain 
I u (r,z)/I u (0,z) = 1 + (l/2)C w (z)r 2 + thus C u {z) 
gives the concavity of the intensity profile at r = 0. 
For example, C u (z) = —A/s^z) for the Gaussian pro- 
file a w (r, z) = exp[-r 2 /s 2 (z)]. 

If our wavepacket E(r, z, t) has the form of enveloped 
oscillations of a certain carrier frequency ujq, then a u (z) 
is a narrow function of frequency around coq (Au/ujo <C 1, 
Au being, e.g., the half-width at 1/e maximum am- 
plitude). The phase and group velocities at which the 
carrier oscillations and the envelope propagate in vac- 
uum are given, respectively, by v p — uJo/d z [uJoz/c + 
<j)oj {z)] and v g = l/d z [ujz/c + 4>u(z)Y Wo , where the prime 
sign means differentiation with respect to w. It is 
then seen that discrepancies of v p and v g from c are 
due, respectively, to axial and frequency dependence 
of Gouy's phase, dependences that originate, also re- 
spectively, from transversal limitation of the wave (i.e., 
from diffraction) and the frequency-dependent nature 
of the diffraction phenomenon. Negative group veloci- 
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ties (with the physical sense discussed by Wang et al. 
||) are in principle allowed by omitting the absolute 
value in the rigorous definition of group velocity jl. 
Since d z (f>(z) — (c/2u})C u> (z), the alternative expressions 
v p /c = [1 + ((? /2u)$)C Uo (z)]~ X , a nd 

7 = {i-^[cu,W-woCi ( 2! )]} , (2) 

can be readily found. The values of v p and v g at a 
given cross section z of a pulsed beam can be then pre- 
dicted from the transversal profile of the monochromatic 
light beam at ujq and neighbouring frequencies. It is 
worthwhile to note that the above expressions of v p and 
v g can be formally obtained from the well-known for- 
mulas v p /c = n~l and v g /c = (n Wo + Uqu'^) -1 for a 
plane pulse in a material medium with refraction index 
n u (z) = l + (c t /2u 2 )C u (z). 

In the most frequent and interesting case of C u (z) < 
(maximum amplitude on- axis), the condition of super- 
luminality [vg(z) > c] is found to be, from Eq. (|^), 
IC^oO 2 )!' > \C uo (z)\/ojq. Negative v g (z) can even take 
place if |C Wo (z)|' > \C LOo (z)\/ujo + 2u;o/c 2 . In terms of the 
equivalent refraction index, superluminality corresponds 
to n' UJg < (1 — n^/ujQ > 0, i.e., to weak normal dis- 
persion or anomalous dispersion, and negative v g corre- 
sponds to n' UJo < —no/ujQ < 0, i.e., to strong anoma- 
lous dispersion. It appears, in conclusion, that the only 
requirement to produce a desired superluminality is to 
introduce a suitable variation with frequency of the con- 
cavity of the transversal profile of the monochromatic 
components. These variations can be easily generated 
in laboratory using simple optical elements as lenses, or 
suitably designed filters: 

Consider an arbitrary pulsed beam E(r, t) at the en- 
trance plane of a lens whose focal length is /, and 
let write its spectrum in the convenient form E u (r) = 
Pujbuj(r), where P u can be identified with the pulse spec- 
trum at a typical point on the lens (e.g., r = 0). The 
dependence of (r) with uj accounts for possible changes 
in the transversal amplitude and phase profile with fre- 
quency. At the focal plane, the spectrum amplitude is, 
as is well-known, a w oc \b u (uir/2Trcf)\, where 

/>oo 

b u (p)=2ir drr6 w (r)J (27rpr) (3) 
Jo 

is the two-dimensional Fourier transform of b^r), Jo(-) 
is the zero-order Bessel function of first class, p = 
(£ 2 + rj 2 ) 1 / 2 , (£,77) being conjugate variables of (x,y). 
The concavity at the focus can be shown to be \C U \ — 
(^/2 7 rc/) 2 |(A x |L|/|L|)p=o|, with A ± = d i6 + d m . Ele- 
mental properties of the Fourier transform lead to the re- 
lationship |(A_l|& w |/|6 w |) p= o| = 47r 2 (r%, where (r 2 ) w = 
[ / °° drr 3 |6^(r)|] / [ f °° drr|b w (r)|] is the squared rms 
width of the transversal profile amplitude on the 

lens. In conclusion, \C U \ = (r 2 ) ^(uj 2 / 'c 2 f 2 ) at the focus. 



In the case of an amplitude profile \b(r)\ independent 
of frequency, the proportionality \C U \ oc uj 2 due to the 
focusing process makes condition |C Wo |' > \C Uo \/u)q of 
superluminality to be satisfied. Indeed, from Eq. (|^), 
the simple result v g /c = (1 — (r 2 ) /2f 2 )~ 1 > I for the 
group velocity at the focus is obtained. The Gaussian 
case b(r) = exp(— r 2 / s 2 ) with plane pulse front and in the 
limit of infinite Fresnel number N Uo = s 2 /\ Ua f (A^ = 
2ttc/ujo) was considered in Ref. [| We see, however, that 
the superluminality at the focus is independent of the 
transversal amplitude and phase profiles, and of A^ . 

Strong superluminality (v g c, v g < 0) may be ob- 
tained if the transversal profile b u {r) on the lens depends 
on frequency. In this case the concavity \C U \ at the fo- 
cus varies with u> due to focusing (oc lu 2 ) and also to the 
possilbe variation of the width (r 2 ) u on the lens with u>. 
From Eq. (j|) the group velocity at the focus is given by 

V 9 = (, 1 <S%0 ^ (OLo V 1 U] 

c \ 2 f 2 2 f 2 J ' W 

which is (as in the previous case) independent of the exact 
shape of b u {r), and also of its possible distortion when 
widening with frequency. Superluminality is enhanced 
with respect to the case of ^-independent profile if (r 2 ) u 
at the lens grows with u> (i.e., {r 2 )' Ug > 0) 

Growth of width with frequency must be obtained 
making use of the spatially dispersive properties of cer- 
tain optical systems. A fair possibility is the use of 
graded mirrors or transparencies Q made of a stack of 
dielectric layers, ]To[ | like the rather simple one shown in 
Fig. [l]a. Such a device is designed to produce a cer- 
tain transversal profile |& WR (r)| of a certain width (r 2 ) Wfl 
when illuminated with a specific frequency lur. When 
illuminated with neighboring frequencies, however, we 
have found that the produced profile changes, and in 
particular its spot size grows significantly with frequency. 
Fig. 0b shows the frequency-dependent width of the out- 
put amplitude profile when the multilayer is illuminated 
with Gaussian beams of different frequencies but con- 
stant width. 

Given uj , (r 2 ) Wo > A 2 , and (r 2 )^ o produced by a 
spatially dispersive system [as those of Fig. 0b in the 
case of the graded mirror] , the focal length to be used to 
obtain v g — 00 at the focus is, from Eq. (0), 

/ 2 = ((r 2 ) Wo +c (r 2 )L )/2 (v g = oo) (5) 

(slightly smaller / will yield a negative group velocity). 
This focal length must moreover verify, for paraxiality, 

/> v/FLo' orfrom ©. 

(rX » (r 2 )^„/uo, (6) 

which can be regarded as the condition for the device 
in order to obtain strong superluminality. In the case 
of the multilayer of Fig. |l|, for the selected uiq = 3.409 
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0.17 mm 2 , and (r 2 V, = 1.645 
J) is loosely satisfied {{r 2 )' UJo ~ 
The focal length for v g = oo is in fact 
0.41 mm. 



/ = 1.7 mm > ^J{r 2 ) u 

Due to the appreciable variation of the width with 
frequency at the lens plane imposed by Eq. (||), we 
must ensure that not only the monochromatic compo- 
nent at w , but at any w within the frequency band 
2Aw, focuses paraxially, i.e., A w <C \J (r 2 )^ <C / for all 
w in 2Aw. If we assume, for simplicity, the approxi- 
mate linear variation of the width J (r 2 ) ~ J (r 2 ) + 

ujq) within 2 Aw, we must then require 

A Wf) for the smallest width 

<C / for the largest 




in the spectrum, and 



These two conditions arc seen to be satisfied if 



Aw < 



(7) 



Eq. (Q) then determines the spectral bandwidth of the 
input pulse to achieve strong superluminality by paraxial 
focusing. For the proposed filter, Eq. (Q) yields Aw <C 
0.206 fs . In practice, the entire spectrum should lie 
within the operation range w > lor of the filter; one can 
then take, e.g., Aw < 0.206/16, which implies an input 
pulse duration At > 155 fs. 

To investigate in more detail the effects associated 
to infinite or negative v g , we have simulated the fo- 
cusing of the pulsed beam E(r, t) whose spectrum 
E u (r) = P u cxp(— r 2 j s\ u) ) has a frequency-dependent 
width sl <uj . Using standard rules of Gaussian beam 
propagation |Q, the spectrum amplitude at any propa- 
gation distance z beyond the lens is given by a^(r, z) = 
[sL^Puj/suiz)} exp[-r 2 /s 2 (z)], where s u {z) = s ijW [(l - 
zj j) 2 + (2cz/ujs 2 Luj ) 2 ] 1 ^ 2 is the Gaussian width at dis- 
tance z, and the spectrum phase is given by <j> u (r,z) = 
-tan" 1 [7riV w (z - f)/z] - tt/2 + u;r 2 /2cR UJo (z), where 
N u = sj^^/X^f is the Gaussian Fresnel number, and 
l/R U0 {z) = [sl^/sUzMI - z/f)/f + (2c/los l ^) 2 z}. 
The group velocity along the z-axis turns out to be, from 
Eq. (0), 
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[1 + (frN u 



(8) 



where a = 1 + 2u} s'r /sl,lu and £, = (z — f)/z. At the 



focus (£ = 0), Eq. (||) reduces to Eq. (0) with the identi- 
fication s 2 L Uq = {r 2 ) UQ . The curve v g (z) is shown in Fig. 
^ for typical sets of parameters. The form of the curve is 
solely determined by N Uo . A maximum is reached at an 
intermediate point between the focus z — f and the waist 
position zq = / /[I + I/t^A^J of the Gaussian beam with 
frequency wo. If N UQ is large enough, the asymmetry of 
the curve v g (z) with respect to z = f becomes inappre- 
ciable. The maximum enhances with increasing conver- 
gence angle sl,uj / /■ Moreover, for given SL,u /f ( as m 



the three curves in Figs. ||a and b), arbitrary enhance- 
ments are achieved by increasing s' T , , . For instance, in 
Fig.|b, ui s' L uo /s LjU)0 = 600 (so that the lhs of H is 2/3 
the rhs). Then an extended region with negative v g ap- 
pears around the focus. This region plays a similar role 
as the gas cell in the experiments of Rcf. S. 

The behavior of a pulse with such an unusual v g is 
more easily understood from Fig. ||a. The solid curve 
shows the arrival time of the pulse peak (that travels 
at velocity v g ) at axial position z, as determined by the 
pulse front equation Q t = [wz/c+0 w (O, z)]' Uo . The slope 
of the curve yields the group velocity in units of c. At 
za and zb, v g is infinite, and negative between them. As 
a consequence (as is apparent from the curve), the pulse 
peak may arrive earlier at axial points behind the focus 
than at points before the focus, e.g., earlier at zb than at 
za- The infinite group velocity at Zb can be associated 
to the fact that the pulse peak appears to arrive instantly 
at zb from zc- Indeed, in the interval (tc,tA), the pulse 
peak arrives at three positions at the same time. 

Of course, the above interpretation for the behavior 
of the pulse peak is valid if the pulse temporal form 
does not experience significant deformation during prop- 
agation, in which case it is also valid for any tempo- 
ral feature of the pulse. In Fig. |3|b we have numer- 
ically verified that this is the case. The pulsed beam 
of spectrum on the lens E u (r) — P^ exp(— r 2 /s| J), 

with P w = ^{b/2)exp[-b 2 (uj - w ) 2 /2] [i.e., P(t) = 
exp(— t 2 /b 2 )] was propagated behind the lens according 
to the Gaussian beam rules (7j, and the time-domain 
field was obtained from numerical integration of Eq. (Q), 
with w = 1.571 fs" 1 , s L , W0 = 0.096 mm, / = 1.92 
mm, and s' T , , = 36.67 mm fs (to obtain, as in Figs, 
gb and |3|a, the values N wo = 4, s L ^ /f = 0.05 and 
^os'l uj / s l,ujq = 600). Condition (^) is satisfied taking 
Aw = 1.309 x 10~ 3 fs" 1 , or b = 2/ Aw = 1528 fs. Fig. ||> 
shows the normalized pulse temporal form at za = 1.861 
mm and zb = 1.974 mm (points A and B of Fig.|]a). 
There is no appreciable deformation in the pulse, which 
arrives at zb (after the focus) 109 fs earlier than at za 
(before the focus). At other distances z, we have also 
observed pulse form invariance. The arrival time of pulse 
peak at each selected z is shown (small squares) in Fig. 
|a. 

Fig. [| offers a more complete off-axis picture of 
the propagation. The pulse front surface t — [luz/c + 
4>u{ r i Z )]L > or sur face formed by all points of equal peak 
time, is depicted at some instants of time. When the 
convergent front advancing towards the focus intersects 
the axis at zc, a new elliptical branch is born at posi- 
tion zb beyond the focus (curves 1), starting its motion 
with v g = oo. At t = 6400 (1 fs later than curve 3), 
the two branches join, transforming into a planar front 
that coincides with the focal plane and backpropagates 
at v g ~ —2c, and a nearly elliptical front around the fo- 
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cus. At longer times (curves 4-6), the situation is nearly 
reversed; in particular, a branch of the front implodes 
at za with v g = oo. In short, the pulse peak (or any 
other pulse feature) diverges from the axial point zg > f 
before it reaches the focal plane. Indeed, the converging 
front never reaches this plane, but dies out at za < /• 

According to our analysis, strong superluminality in 
vacuum, including negative group velocities, could be ex- 
perimentally demonstrated by focusing rather arbitrary 
pulses, under the only condition that the spot size varies 
with frequency appropriately. The required variation 
has been shown to be loosely attained with a multilayer 
graded mirror. Other devices are under study. Superlu- 
minality is here an effect of the dephasing between the 
different spectral components of the pulse due to the fre- 
quency dependence of Gouy's phase. 

The authors thank A. Piegari for helpful discussions. 
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FIG. 1: a) An all-dielectric multilayer high transmittance 
coating, obtained from a stack of quarter -wave (Ahj^ — 555 
nm, or ljr = 3.396 fs -1 ) layers of alternate high (jib = 
2.20) and low index (jll = 1.48) on a transparent substrate 
(ns = 1-48), producing a radially variable transmittance. b) 
The radius of the multilayer is r = 10 mm. We used the 
simulation computer program of Ref. to obtain the ra- 
dial transmittance curve t^(r) for ujr and higher frequencies. 
When the multilayer is illuminated with the frequency in- 
dependent Gaussian profile exp(— r 2 /so), so = 0.5 mm, the 
output is t^(r) exp(— r 2 /sg), whose width is depicted (open 
circles). For convenience we choose ujo = 3.409 fs -1 , for which 
{r 2 )^ = 0.17 mm 2 and (r 2 )i = 1.645 mm 2 fs 

FIG. 2: Group velocity of focused Gaussian pulsed beams as 
functions of propagation distance from the lens. In all cases 
N Uo = 4, s LM0 /f = 0.05. (a) u s' L ^Js L ^ = 0,0.003, for 
the solid and dashed curves respectively, (b) uJos' L uJa /sL^ = 
600. 

FIG. 3: (a) Position of the pulse peak as it propagates, as pre- 
dicted by t = [ujz/c + <^(0, z)]L and numerically calculated, 
(b) Normalized pulse temporal form at propagation distances 
given. See the text for numerical values of the parameters. 

FIG. 4: Pulse front at times a) 1: 6346 fs, 2: 6378 fs, 3: 6399 
fs. b) 4: 6401 fs, 5: 6422 fs, 6: 6454 fs. 
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